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ABSTRACT 

Recently, the High Resolution Fly's Eye (HiRes) experiment claims that there is no small scale anisotropy in the 
arrival distribution of ultra-high energy cosmic rays (UHECRs) above E > 10 19 eV contrary to the Akeno Giant 
Air Shower Array (AGASA) observation. In this paper, we discuss the statistical significance of this discrepancy 
between the two experiments. We calculate arrival distribution of UHECRs above 10 19 eV predicted by the source 
models constructed using the Optical Redshift Survey galaxy sample. We apply the new method developed by 
us for calculating arrival distribution in the presence of the galactic magnetic field. The great advantage of this 
method is that it enables us to calculate UHECR arrival distribution with lower energy (~ 10 19 eV) than previous 
studies within reasonable time by following only the trajectories of UHECRs actually reaching the earth. It has 
been realized that the small scale anisotropy observed by the AGASA can be explained with the source number 
density ~ 10" 5 ~" 6 Mpc" 3 assuming weak extragalactic magnetic field (B < 1 nG). We find that the predicted 
small scale anisotropy for this source number density is also consistent with the current HiRes data. We thus 
conclude that the statement by the HiRes experiment that they do not find small scale anisotropy in UHECR 
arrival distribution is not statistically significant at present. We also show future prospect of determining the 
source number density with increasing amount of observed data. 

Subject headings: cosmic rays — methods: numerical — ISM: magnetic fields — galaxies: general — 
large-scale structure of universe 



1. INTRODUCTION 

The small scale anisotropy in the observed arrival distribu- 
tion of ultra-high energy cosmic rays (UHECRs) is key in iden- 
tifying their sources which will tell us a great deal about accel- 
eration mechanisms, composition of UHECRs, and so on. The 
Akeno Giant Air Shower Array (AGASA) observation reveals 
the existence of the small scale clusterings in the isotropic ar- 
rival distribution of UHECRs (Takeda et al. 1999, 2001). The 
current AGASA data set of 57 events above 4 x 10 19 eV con- 
tains four doublets and one triplet within a separation angle 
of 2.5°. Chance probability to observe such clusters under an 
isotropic distribution is only about 1 % (Hayashida et al. 2000; 
Takeda et al. 2001). 

On the other hand, the recent claim by the High Resolution 
Fly's Eye (HiRes; Wilkinson et al. 1999) experiment is that 
there is no small scale anisotropy in the UHECR arrival distri- 
bution above E > 10 19 eV (162 events) which is observed by 
a stereo air fluorescence detector (Finley et al. 2003). The dis- 
crepancy between the AGASA and HiRes involves a number 
of consideration. Among them are the different numbers of the 
events, the possibility of the measured energies being shifted 
between the two experiments, (Marco, Blasi, & Olinto 2003) 
and the different angular resolutions. The purpose of this paper 
is to discuss the first possibility. That is, we consider whether 
the discrepancy between the two experiments is statistically sig- 
nificant or not at present. 

There is furthermore the contradiction on the presence or ab- 
sence of the GZK cutoff (Greisen 1966; Zatsepin & Kuz'min 
1966) in the cosmic-ray energy spectrum due to photopion pro- 
duction with the photons of the cosmic microwave background 
(CMB). The HiRes spectrum shows the GZK cutoff (Abu- 



Zayyad et al. 2002), but the AGASA spectrum does not (Takeda 
et al. 1998). The statistical significance of this discrepancy is 
already discussed in Marco, Blasi, & Olinto (2003) and shown 
to be low at present. This issue is left for future investigation by 
new large -aperture detectors under development, such as South 
and North Auger project (Capelle et al. 1998), the EUSO (Ben- 
son & Linsley 1982), and the OWL (Cline & Stecker 2000) 
experiments. 

When we consider the statistical significance of the small 
scale anisotropy, we need to compare the observations with 
numerical calculations predicted by some kind of source dis- 
tributions. In this paper, we use the Optical Redshift Survey 
(ORS; Santiago et al. 1995) galaxy sample to construct realis- 
tic source models of UHECRs. This sample has been adopted in 
our previous studies (Yoshiguchi et al. 2003a; Yoshiguchi, Na- 
gataki, & Sato 2003c,d). It has been realized that the small scale 
anisotropy well reflect the number density of UHECR sources 
(Yoshiguchi et al. 2003a; Yoshiguchi, Nagataki, & Sato 2003c; 
Blasi & Marco 2004). And also, it is shown that the small scale 
anisotropy above E = 4 x 10 19 eV observed by the AGASA 
can be explained with the source number density ~ io _5 ~- 6 
assuming weak extragalactic magnetic field (B < 1 nG, EGMF) 
and that observed UHE particles are protons (Yoshiguchi et al. 
2003a; Yoshiguchi, Nagataki, & Sato 2003c; Blasi & Marco 
2004). We thus take the source number density as a parameter 
of our source model, and discuss the statistical significance of 
the small scale anisotropy by considering to what extent we can 
determine the source number density from the comparison of 
our model predictions with the HiRes observation. 

As mentioned above, the claim by the HiRes experiment on 
the small scale anisotropy is based on the data of UHECR ar- 
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rival directions above 10 19 eV. UHE protons at such energies 
are significantly deflected by the galactic magnetic field (GMF) 
by a few - ~ 10 degrees. In order to accurately calculate the 
expected UHECR arrival distribution and compare with the ob- 
servations, this effect should be taken into account. 

In our recent work (Yoshiguchi, Nagataki, & Sato 2003d), 
we presented a new method for calculating UHECR arrival dis- 
tribution which can be applied to several source location sce- 
narios including modifications by the GMF. Brief explanation 
of our method is as follows. We numerically calculate the prop- 
agation of anti-protons from the earth toward the outside of the 
Galaxy (we set a sphere centered around the Galactic center 
with radius r src =40 kpc as the boundary condition), consider- 
ing the deflections due to the GMF. The anti-protons are ejected 
isotropically from the earth. By this calculation, we can obtain 
the trajectories and the sky map of position of anti-protons that 
have reached to the boundary at radius r src = 40 kpc. 

Next, we regard the obtained trajectories as the ones of PRO- 
TONS from the outside of the galaxy toward the earth. Also, 
we regard the obtained sky map of position of anti -proton at the 
boundary as relative probability distribution (per steradian) for 
PROTONs to be able to reach to the earth for the case in which 
the flux of the UHE protons from the extra-galactic region is 
isotropic (in this study, this flux corresponds to the one at the 
boundary r src = 40 kpc). The validity of this treatment is sup- 
ported by the Liouville's theorem. When the flux of the UHE 
protons at the boundary is anisotropic (e.g., the source distribu- 
tion is not isotropic), this effect can be included by multiplying 
this effect (that is, by multiplying the probability density of ar- 
rival direction of UHE protons from the extra-galactic region at 
the boundary) to the obtained relative probability density distri- 
bution mentioned above. 

By adopting this new method, we can consider only the tra- 
jectories of protons which arrive to the earth, which, of course, 
helps us to save the CPU time efficiently and makes calcula- 
tion of propagation of CRs even with low energies (~ 10 19 eV) 
possible within a reasonable time. 

With this method, we calculate the arrival distribution of 
UHE protons above 10 19 eV for our source models. We also 
consider the energy loss processes when UHE protons propa- 
gate in the intergalactic space. Using the two point correlation 
function as statistical quantities for the small scale anisotropy, 
we compare our model prediction with the HiRes observation. 
We then discuss to what extent we can determine the source 
number density of UHECRs by the current HiRes data. We 
also briefly discuss future prospect of determining the source 
number density with the event number expected by future ex- 
periments. 

In section 2, we introduce the GMF model. We explain the 
method for calculating UHECR arrival distribution in section 3. 
Results are shown in section 4. In section 5, we summarize the 
main results. 

2. GALACTIC MAGNETIC FIELD 

In this study, we adopt the GMF model used in Alvarez- 
Muniz, Engel & Stanev (2002); Yoshiguchi, Nagataki, & Sato 
(2003d, 2004), which is composed of the spiral and the dipole 
field. In the following, we briefly introduce this GMF model. 

Faraday rotation measurements indicate that the GMF in the 
disk of the Galaxy has a spiral structure with field reversals at 
the optical Galactic arms (Beck 2001). We use a bisymmetric 
spiral field (BSS) model, which is favored from recent work 



(Han, Manchester, & Qiao 1999; Han 2001). The Solar Sys- 
tem is located at a distance r\\ =R® = 8.5 kpc from the center 
of the Galaxy in the Galactic plane. The local regular mag- 
netic field in the vicinity of the Solar System is assumed to be 
BsoVdi ~ 1.5 /iG in the direction / = 90°+ p where the pitch angle 
is p = -10° (Han & Qiao 1994). 

In the polar coordinates {r\ \ , 4>), the strength of the spiral field 
in the Galactic plane is given by 



B(r\\,(f>)=B 



^ ) cos[ 0-/3k3 
r \\ ) V r o 



(1) 



where B Q = 4.4 /j,G, r Q = 10.55 kpc and (3 = l/tmp = -5.67. 
The field decreases with Galactocentric distance as 1 jr\ and it 
is zero for n i > 20 kpc. In the region around the Galactic center 
(ni < 4 kpc) the field is highly uncertain, and thus assumed to 
be constant and equal to its value at r\ \ =4 kpc. 

The spiral field strengths above and below the Galactic plane 
are taken to decrease exponentially with two scale heights 
(Stanev 1996) 

id/- a m lor ,am / ex P(" z ) : |z|<0.5kpc 
l*(ni^)l = |B(ni^)l( exp(^)exp(f): \z\ > 0.5 kpc 

(2) 

where the factor exp(-3/8) makes the field continuous in z- The 
BSS spiral field we use is of even parity, that is, the field direc- 
tion is preserved at disk crossing. 

Observations show that the field in the Galactic halo is much 
weaker than that in the disk. In this work we assume that the 
regular field corresponds to a AO dipole field as suggested in 
(Han 2002). In spherical coordinates (r,6,(p) the (x,y,z) com- 
ponents of the halo field are given by: 

B x = —2> /iG sin # cos # cos (^/V 3 



By =-3 sin 9 cos 6 sin (p/r 3 
B z = fi G (1- 3 cos 2 9)/r 3 



(3) 



where /ig ~ 184.2 /iG kpc 3 is the magnetic moment of the 
Galactic dipole. The dipole field is very strong in the central 
region of the Galaxy, but is only 0.3 /iG in the vicinity of the 
Solar system, directed toward the North Galactic Pole. 

There may be a significant turbulent component, B mn , of 
the Galactic magnetic field. Its field strength is difficult to 
measure and results found in literature are in the range of 
B mn = 0.5 . . .2B reg (Beck 2001). However, we neglect the ran- 
dom field throughout the paper for simplicity. Possible depen- 
dence of the results on the assumption is discussed in the final 
section. 

3. NUMERICAL METHOD 
3.1. Propagation of UHE protons in the Intergalactic Space 

The energy spectrum of UHECRs injected at extragalactic 
sources is modified by the energy loss processes when they 
propagate in the intergalactic space. In this subsection, we ex- 
plain the method of Monte Carlo simulations for propagation 
of UHE protons in intergalactic space. 

We assume that the composition of UHECRs are protons 
which are injected with a power law spectrum within the range 
of (10 19 - 10 22 )eV. We inject 10000 protons in each of 31 en- 
ergy bins (10 bins per decade of energy). Then, UHE protons 
are propagated including the energy loss processes (explained 
below) over 3 Gpc for 15 Gyr. We take a power law index as 
2.6 to fit the predicted energy spectrum to the one observed by 
the HiRes (Marco, Blasi, & Olinto 2003). 
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UHE protons below ~ 8 x 10 eV lose their energies mainly 
by pair creations and adiabatic energy losses, and above it by 
photopion production (Berezinsky & Grigorieva 1988; Yoshida 
& Teshima 1993) in collisions with photons of the CMB. We 
treat the adiabatic loss as a continuous loss process. We cal- 
culate the redshift z of source at a given distance using the 
cosmological parameters H = 71 km s" 1 Mpc" 1 , Q m = 0.27, 
and Oa = 0.73. Similarly, the pair production can be treated 
as a continuous loss process considering its small inelasticity 
(~ 10~ 3 ). We adopt the analytical fit functions given by Chodor- 
owski, Zdziarske, & Sikora (1992) to calculate the energy loss 
rate for the pair production on isotropic photons. 

On the other hand, protons lose a large fraction of their en- 
ergy in the photopion production. For this reason, its treatment 
is very important. We use the interaction length and the en- 
ergy distribution of final protons as a function of initial proton 
energy which is calculated by simulating the photopion produc- 
tion with the event generator SOPHIA (Mucke et al. 2000). The 
same approach has been adopted in our previous studies (Ide et 
al. 2001; Yoshiguchi et al. 2003a; Yoshiguchi, Nagataki, & Sato 
2003c,d). 

In this study, we neglect the effect of the EGMF because 
of the following two reasons. First, numerical simulations 
of UHECR propagation in the EGMF including lower energy 
(~ 10 19 eV) ones take a long CPU time. Secondly, we show 
in our previous study that small scale clustering observed by 
the AGAS A can be well reproduced in the case of weak EGMF 
(B < InG) (Yoshiguchi et al. 2003a; Yoshiguchi, Nagataki, & 
Sato 2003c). The purpose of this paper is to discuss whether 
the claim by the HiRes experiment is contrary to the AGASA 
result with statistical significance. We should calculate the ar- 
rival distribution of UHECRs in the situation that the AGASA 
experiment can be explained. 

Of course, effects of strong EGMF (say, ^100 nG) on 
UHECR propagation are intensively investigated (Sigl, Miniati, 
& Ensslin 2003,b, 2004; Aloisio & Berezinsky 2004), and can 
be hardly excluded considering the limited amount of observa- 
tional data on the EGMF. However, we assume extremely weak 
EGMF throughout the paper. 

3.2. Source Distribution 

In this study, we apply the method for calculating the 
UHECR arrival distribution with modifications by the GMF 
(section 3.3) to our source location scenario, which is con- 
structed by using the ORS (Santiago et al. 1995) galaxy catalog. 

In order to calculate the distribution of arrival directions of 
UHECRs realistically, there are two key elements of the galaxy 
sample to be corrected. First, galaxies in a given magnitude- 
limited sample are biased tracers of matter distribution because 
of the flux limit (Yoshiguchi et al. 2003b). Although the sample 
of galaxies more luminous than -20.5 mag is complete within 
80 h~ x Mpc (where h is the Hubble constant divided by 100 km 
s" 1 and we use h = 0.71), it does not contain galaxies outside it 
for the reason of the selection effect. We distribute sources of 
UHECRs outside 80 h~ l Mpc homogeneously. Their number 
density is set to be equal to that inside 80 h~ x Mpc. 

Secondly, our ORS sample does not include galaxies in the 
zone of avoidance (\b\ < 20°). In the same way, we distribute 
UHECR sources in this region homogeneously, and calculate 
its number density from the number of galaxies in the observed 
region. 

As mentioned above, we take the number density of UHECR 



sources as our model parameter when we construct the source 
models. For a given number density, we randomly select galax- 
ies from the above sample with probability proportional to the 
absolute luminosity of each galaxy. In this paper, the model 
parameters are taken to be 10" 3 , 10" 4 , 10" 5 , and 10" 6 Mpc" 3 . 

3.3. Calculation of the UHECR Arrival Distribution with 
modifications by the GMF 

This subsection provides the method of calculation of 
UHECR arrival distribution with modifications by the GMF. 
The explanations largely follows our recent paper (Yoshiguchi, 
Nagataki, & Sato 2003d). We start by injecting anti-protons 
from the earth isotropically, and follow each trajectory until 

1 . anti -proton reaches a sphere of radius 40 kpc centered at 
the galactic center, or 

2. the total path length traveled by anti-proton is larger than 
200 kpc, 

by integrating the equations of motion in the magnetic field. 
It is noted that we regard these anti-protons as PROTONs in- 
jected from the outside of the Galaxy toward the earth. The 
number of propagated anti-proton is 2,000,000. We have 
checked that the number of trajectories which are stopped by 
the limit (2) is smaller than 0.1% of the total number. The en- 
ergy loss of protons can be neglected for these distances. Ac- 
cordingly, we inject anti-protons with injection spectrum simi- 
lar to the observed one ~ E~ 2 1 . (Note that this is not the energy 
spectrum injected at extragalactic sources.) 

The result of the velocity directions of anti-protons at the 
sphere of radius 40 kpc is shown in the right panel of figure 1 
in the galactic coordinate. From Liouville's theorem, if the 
cosmic-ray flux outside the Galaxy is isotropic, one expects an 
isotropic flux at the earth even in the presence of the GMF. This 
theorem is confirmed by numerical calculations shown in figure 
6 of Alvarez-Muniz, Engel & Stanev (2002), which is the same 
figure as our figure 1 except for threshold energy. Thus, the 
mapping of the velocity directions in the right panel of figure 1 
corresponds to the sources which actually give rise to the flux 
at the earth in the case that the sources (including ones which 
do not actually give rise to the flux at the earth) are distributed 
uniformly. 

We calculate the UHECR arrival distribution for our source 
scenario using the numerical data of the propagation of UHE 
anti-protons in the Galaxy. Detailed method is as follows. At 
first, we divide the sky into a number of bins with the same 
solid angle. The number of bins is taken to be 360(0 x 200(b). 
We then distribute all the trajectories into each bin according 
to their directions of velocities (source directions) at the sphere 
of radius 40 kpc. Finally, we randomly select trajectories from 
each bin with probability Pse\ ec (j ,k,E) defined as 

PselecO,A:,£)oC^^ ^=2/7 ' ^ 

Here subscripts j and k distinguish each cell of the sky, di is 
distance of each galaxy within the cell of (j,k), and the summa- 
tion runs over all of the galaxies within it. E is the energy of 
proton, and dN/dE(di,E) is the energy spectrum of protons at 
our galaxy injected at a source of distance di. 

The normalization of P se \ ec (j,k,E) is determined so as to set 
the total number of events equal to a given number, for example, 
the event number of the current HiRes data. When P se i e c > 1, 
we newly generate events with number of (P se iec - 1) x N(j,k), 
where N(j,k) is the number of trajectories within the sky cell 
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FIG. 1. — Arrival directions of anti-protons with E > 10 eV at the sphere of Galactocentric radius of 40 kpc (right panel) in the galactic coordinate. The 
anti-protons are injected at the earth isotropically (as shown in the left panel) with an injection spectrum E~ ln . 



of ( j,k). The arrival angle of newly generated proton (equiva- 
lently, injection angle of anti-proton) at the earth is calculated 
by adding a normally distributed deviate with zero mean and 
variance equal to the experimental resolution 2.8° (1.8°) for 
E > 10 19 eV (4 x 10 19 eV) to the original arrival angle. 

3.4. Statistical Methods 

A standard tool for searching the small scale anisotropy is 
the two point correlation function. This subsection provides the 
explanation of this function. 

We start from a set of simulated events. For each event, we 
divide the sphere into concentric bins of angular size AO, and 
count the number of events falling into each bin. We then divide 
it by the solid angle of the corresponding bin, that is, 

2tt cos6>-cos(# + A60 ' 

where 4> denotes the separation angle of the two events. AO is 
taken to be 2° in this analysis. We also evaluate N un i(9) for an 
equal number of events simulated in a Monte Carlo for uniform 
source distribution. The estimator for the correlation function 
is then defined as w(0) = N(0)/N uni (0)-l. 

The two point correlation function of the HiRes data with 
164 events above 10 19 eV do not show significant small scale 
anisotropy, and is consistent with isotropic source distribution 
within la confidence level. We thus compare the two point 
correlation functions predicted by our source models with that 
expected for isotropic source distribution rather than the HiRes 
data itself. When we consider to what extent we can determine 
the source number density of UHECRs, we have to quantify 
deviations of predictions by our source models from isotropic 
sources. To do so, we define the variable \ 2 as 

2 = J_ (Wthffl)-Wuniffi)) 2 

X AW (a\ h ) 2 + (ar) 2 ' 

where 8\ is the angle at i-th bin. of 1 and cr" nl are the statistical 
fluctuations of w t h and w un \ at the angle 0\. The summation is 
taken over the angular bins, whose number is denoted by AW 
Even if we specify the source number density, source dis- 
tribution itself can not be determined because of randomness 



when we select galaxies from the ORS sample. Therefore, we 
first evaluate x 2 f° r a source distribution in the case of a given 
number density of UHECR sources. We then repeat such cal- 
culations for a number of realizations of the source selection, 
and obtain the probability distribution of \ 2 f° r a source num- 
ber density. The source selections are performed 100 times for 
all the source number densities. 

4. RESULTS 

4. 1 . Arrival Distribution of UHECRs 

In this subsection, we present the results of arrival distribu- 
tion of UHECRs above 10 19 eV obtained by using the method 
explained in section 3.3. Figure 2 and 3 show realizations of 
the event generations in the case of the source number density 
10~ 3 and 10~ 5 Mpc" 3 , respectively. The events are shown by 
color according to their energies. Note that the event number of 
HiRes data is 162 (E > 10 19 eV). This roughly corresponds to 
300 events of figure 2 and 3, considering the difference of the 
range <5 (declination) between the observation and numerical 
calculation. 

As is evident from figure 2 and 3, arrival distributions of 
UHECRs are quite isotropic. The harmonic analysis to the right 
ascension distribution of events is the conventional method to 
search for global anisotropy of cosmic ray arrival distribution. 
Using this analysis, we have checked that large scale isotropy 
predicted by our source models is consistent with isotropic 
source distribution within a 90% confidence level with event 
number equal to the current HiRes and AGASA observations. 
On the other hand, arrival distributions observed by the HiRes 
and AGASA do not also show any significant global anisotropy. 
There is no discrepancy between the two experiments. Thus we 
do not discuss large scale isotropy in the following. 

Comparing figure 2 and 3, we easily find that arrival dis- 
tribution expected for smaller source number density exhibits 
stronger anisotropy on small angle scale. This is simply be- 
cause contributions from each source to the total events are 
larger for smaller source number density. Furthermore, the 
clustered events are aligned in the sky according to the order of 
their energies, reflecting the direction of the GMF at each direc- 
tion. This interesting feature becomes evident with increasing 
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FIG. 2. — Arrival directions of protons with E > lO 1 ' eV predicted by a source distribution for the source number density 10 3 Mpc 3 in Galactic coordinates. 
Events are shown by color according to their energies. 




FIG. 3. — Same as figure 2, but for the source number density 10 5 Mpc 
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amount of the event number. 



4.2. Statistical significance of the small scale anisotropy 

Next, we discuss the statistics on arrival distribution of UHE- 
CRs. Figure 4 and 5 show the two point correlation function 
expected in the case of source distributions of figure 2 and 3, re- 
spectively. We calculate the two point correlation function for 
the simulated events within only 0° < 5 < 90° in accordance 
with the HiRes exposure. For this reason, the event numbers 
written in figure 4 and 5 differ from the ones in figure 2 and 
3. The error bars represent 1 a statistical fluctuations of our 
numerical calculations due to the finite number of events. The 
shaded regions represent 1 a fluctuations for uniform source 
distribution. The event selections are performed 1000, 100, and 
30 times for the event number 150, 1500, and 5000, respec- 
tively. It is noted that we compare our numerical results with 
not the HiRes data itself but the prediction of uniform source 
distribution. We do not have to consider the non-uniform ob- 
servation in the angular sky. 
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FIG. 4. — Two point correlation function predicted by a realization of the 
source selection for the number density 10~ 3 Mpc -3 . We calculate the two 
point correlation function for the simulated events within only 0° < <5 < 90° 
in accordance with the HiRes exposure. The error bars represent 1 a statistical 
fluctuations due to the finite number of events for our numerical calculations. 
The shaded regions represent 1 <r fluctuations for uniform source distribution. 
The event selections are performed 1000, 100, and 30 times for the event num- 
ber 150, 1500, and 5000, respectively. 

As mentioned before, it is evident that smaller source number 
density predicts stronger correlation between events at small 
angle scale. The small scale anisotropy well reflects the num- 
ber density of UHECR sources. It is already shown that the 
small scale anisotropy above E = 4 x 10 19 eV observed by 
the AGASA can be explained with the source number den- 
sity ~ 10~ 5 6 Mpc ~ 3 assuming weak extragalactic magnetic 
field (B < 1 nG, EGMF) (Yoshiguchi et al. 2003a; Yoshiguchi, 
Nagataki, & Sato 2003c; Blasi & Marco 2004). We thus dis- 
cuss the statistical significance of the discrepancy between the 
HiRes and AGASA by considering to what extent we can de- 
termine the UHECR source number density from the current 
HiRes data. 



Even if we specify the source number density, source dis- 
tribution itself can not be determined because of randomness 
when we select galaxies from the ORS sample. Therefore, we 
first evaluate the statistical significance of deviation of the two 
point correlation function predicted by a source distribution for 
a given number density from that expected for uniform source 
distribution. We then repeat such calculations for a number of 
realizations of the source selection, and obtain the probability 
distribution of x 2 for a given source number density. 
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FIG. 5. — Same as figure 4, but for the number density 10~ 5 Mpc -3 . 
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FIG. 6. — Probability distribution of statistical significance of deviation of 
our model prediction from that of isotropic sources (x 2 ) when the source se- 
lection from our ORS sample is performed 100 times. The data of the two 
point correlation at the smallest angle are used (N bin = 1). 

The results of the probability distribution of x 2 f° r M>in = 1 
and 5 are shown as histograms in figure 6 and 7, respectively. 
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The source selection is performed 100 times for all the source 
number densities. As we can see from the figures, differences 
between probability distributions of y 2 f° r different source 
number densities are clearer when AW = 1 ■ This is related to 
the fact that the two point correlation function is most sensitive 
to the source number density at the smallest angular bin (see, 
figure 4 and 5). We thus focus our attention to the result for 
AW = 1 in the following. 




FIG. 7. — Same as figure 6, but for Afbin = 5. 

As mentioned, we know that the small scale anisotropy ob- 
served by the AGASA can be explained with the number den- 
sity ~ 10" 6 ~" 5 Mpc" 3 . From figure 6, it is found that the small 
scale anisotropy predicted by this source number density is also 
consistent with the prediction of isotropic source distribution 
when the event number is equal to the HiRes data (~ 150). It is 
noted that the AGASA data have about 1000 events above 10 19 
eV. The possibility that the number density is about ~ 10" 6 ~" 5 
Mpc" 3 can not be ruled out, although the HiRes result seems to 
be in agreement with the isotropic source distribution. We thus 
conclude that the statement by the HiRes experiment that they 
do not find small scale anisotropy in UHECR arrival distribu- 
tion is not statistically significant at present. 

We now turn to the future prospects of determining the 
UHECR source number density. The probability distributions 
of x 2 with larger event number than the current observations are 
also shown in figure 6. Of course, there is no observational data 
to be compared. We thus compare the model predictions with 
the case of isotropic source distribution. The results should be 
interpreted as quantifying deviations of our calculations from 
the prediction of isotropic sources in the unit of statistical fluc- 
tuations. 

As is evident from figure 6, difference between the probabil- 
ity distribution of \ 2 f° r different source number densities be- 
come clear when the number of data increases. We will be able 
to determine the source number density by the future observa- 
tions of the small scale anisotropy. In order to discuss to what 
accuracy we can know the source number density, we need to 
calculate the probability distribution of x 2 more precisely with 
a large number of simulations of the source selection. This is- 



sue is beyond the scope of the present paper. We will study it 
in a future publication. 

5. SUMMARY AND DISCUSSION 

In this paper, we considered the statistical significance of the 
discrepancy between the HiRes and the AGASA experiment 
on the small scale anisotropy in arrival distribution of UHE- 
CRs. We calculated arrival distribution of UHECRs above 10 19 
eV predicted by our source models constructed using the ORS 
galaxy sample. We applied the new method developed by us 
(Yoshiguchi, Nagataki, & Sato 2003d) for calculating arrival 
distribution with the modifications by the GMF The great ad- 
vantage of this method is that it enables us to calculate UHECR 
arrival distribution with lower energy (~ 10 19 eV) than previ- 
ous studies by following only the trajectories actually reaching 
the earth. 

It has been realized that the small scale anisotropy in the 
UHECR arrival distribution well reflects the number density 
of UHECR sources (Yoshiguchi et al. 2003a; Yoshiguchi, Na- 
gataki, & Sato 2003c; Blasi & Marco 2004). And also, it is 
shown that the small scale anisotropy above E = 4 x 10 19 eV ob- 
served by the AGASA can be explained with the source number 
density ~ 10~ 5 ~" 6 assuming weak extragalactic magnetic field 
(B < 1 nG, EGMF) and that observed UHE particles are protons 
(Yoshiguchi et al. 2003a; Yoshiguchi, Nagataki, & Sato 2003c; 
Blasi & Marco 2004). We thus took the source number density 
as a parameter of our source model, and discussed the statisti- 
cal significance of the small scale anisotropy by considering to 
what extent we can determine the source number density from 
the comparison of the model predictions with the HiRes obser- 
vation. 

When we consider to what extent we can determine the 
source number density of UHECRs, we have to quantify devia- 
tions of predictions by our source models from the observation. 
And also, even if we specify the source number density of our 
model, source distribution itself can not be determined because 
of randomness when we select galaxies from the ORS sample. 
Therefore, we first evaluated the statistical significance of de- 
viation (x 2 ) of the two point correlation function predicted by 
a source distribution for a given number density from that ex- 
pected for uniform source distribution. We then repeated such 
calculations for a number of realizations of the source selection. 
This gave the probability distributions of \ 2 f° r various number 
densities of UHECR sources. 

The results for AW = 1 and 5 show that differences between 
probability distributions of \ 2 f° r different source number den- 
sities are clearer when AW = L This is related to the fact 
that the two point correlation function is most sensitive to the 
source number density at the smallest angular bin. We thus 
focused our attention to the result of AW = 1 • As mentioned 
above, we know that the small scale anisotropy observed by the 
AGASA can be explained with the number density ~ 10" 6 ~" 5 
Mpc" 3 . It is found that the small scale anisotropy predicted by 
this source number density is also consistent with the predic- 
tion of isotropic source distribution when the event number is 
equal to the HiRes data (~ 150). Note that the data number of 
AGASA is about 1000 above 10 19 eV. The possibility that the 
number density is about ~ 10" 6 ~~ 5 Mpc" 3 can not be ruled out, 
although the HiRes result seems to be in agreement with the 
isotropic source distribution. We thus concluded that the state- 
ment by the HiRes experiment that they do not find small scale 
anisotropy in UHECR arrival distribution is not statistically sig- 
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nificant at present. 

We also discussed the future prospects of determining the 
UHECR source number density. As is evident from figure 6, 
difference between the probability distribution of y 2 for differ- 
ent source number densities become clear when the number of 
data increases. We will be able to determine the source number 
density by the future observations of the small scale anisotropy. 
In order to discuss to what accuracy we can know the source 
number density, we need to calculate the probability distribu- 
tion of x 2 more precisely with a large number of simulations of 
the source selection. This issue deserves further investigation. 

Finally, we mention the assumptions made in this paper. We 
neglected the effects of the extragalactic magnetic field and the 
random component of the GMF. If we take these effect into ac- 
count, the small scale anisotropy obtained by numerical calcu- 
lations become less obvious. In this case, our statement that the 
small scale anisotropy for <~ 10" 5 ~" 6 Mpc" 3 is also consistent 
with the prediction of isotropic sources when the event number 



is equal to the HiRes data (~ 150) remains to be valid. Hence, 
the assumptions are not so important for our conclusion. 

Note added: While we are finishing this study, a paper 
by the HiRes collaboration (Abbasi et al. 2004) appears on 
the net, where they present the results of a search for small 
scale anisotropy in the observed arrival distribution of UHE- 
CRs. They find no small scale anisotropy in 271 events above 
10 19 eV, which is slightly larger than the event number consid- 
ered in the present study (164). However, this increase of the 
event number will not affect our conclusion very much, as seen 
from figure 6. 
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